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ABSTRACT 

We present the first statistical survey of the properties of the 12 CO(1-0) and 
12 CO(3-2) line emission from the nuclei of a nearly complete subsample of 60 
infrared (IR) luminous galaxies selected from SCUBA Local Universe Galaxy 
Survey (SLUGS). This subsample is flux limited at > 5.24 Jy with far-IR 

(FIR) luminosities mostly at Lfir > 10 10 L Q . We compare the emission line 
strengths of 12 CO(1-0) and (3-2) transitions at a common resolution of 15". The 
measured 12 CO(3-2) to (1-0) line intensity ratios r 31 vary from 0.22 to 1.72 with a 
mean value of 0.66 for the sources observed, indicating a large spread of the degree 
of excitation of CO in the sample. These CO data, together with a wide range of 
data at different wavelengths obtained from the literature, allow us to study the 
relationship between the CO excitation conditions and the physical properties 
of gas/dust and star formation in the central regions of galaxies. Our analysis 
shows that there is a non-linear relation between CO and FIR luminosities, such 
that their ratio Lco/Lfir decreases linearly with increasing Lpm- This behavior 
was found to be consistent with the Schmidt Law relating star formation rate 
to molecular gas content, with an index N = 1.4 ± 0.3. We also find a possible 
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dependence of the degree of CO gas excitation on the efficiency of star forming 
activity. Using the large velocity gradient (LVG) approximation to model the 
observed data, we investigate the CO-to-H 2 conversion factor X for the SLUGS 
sample. The results show that the mean value of X for the SLUGS sample is 
lower by a factor of 10 compared to the conventional value derived for the Galaxy 
if we assume the abundance of CO relative to H 2 , Zco = 1CT 4 . For a subset of 
12 galaxies with H I maps, we derive a mean total face-on surface density of 
H 2 +H I of about 42 M Q pc~ 2 within about 2 kpc of the nucleus. This value is 
intermediate between that in galaxies like our own and those with strong star 
formation. 

Subject headings: radio lines: galaxies — galaxies: ISM — galaxies: starburst - 
ISM: evolution — ISM: molecules — surveys 

1. Introduction 

Knowledge of the properties and evolution of the gas and dust content in IR luminous 
galaxies is essential for understanding the cause and temporal evolution of starburst activ- 
ity. In particular, studies of such galaxies in the nearby universe are essential as a step 
in understanding the role of the starburst phenomenon in the cosmic evolution of galaxies. 
Therefore, there is a need to investigate large statistical samples of IR luminous galaxies 
using a multitude of different types of data, including CO, H I, and continuum in the sub- 
millimeter (sub-mm), FIR and radio, in order to constrain theories of how the interstellar 
medium (ISM) evolves (e.g. Eales & Edmunds 1996, 1997). A significant step was achieved 
by Dunne et al. (2000) for the nearby universe by investigating the sub-mm properties of 
the dust in a complete sample of 104 galaxies. The SLUGS survey is based on the Revised 
Bright Galaxy Sample of IRAS galaxies (Soifer et al. 1987) within -10° < 5 < +50° and 
with cz > 1900 km s~ x and a flux limit of S^Ofim > 5.24 Jy. By virtue of the selection, the 
SLUGS galaxies have the FIR luminosity L FIR > 10 10 L , and hence have greater than the 
average star formation rate. A study of molecular gas of a complete flux-limited subsample 
is important to complement the study of the dust. 

Low lying CO rotational line transitions at mm and sub-mm wavelengths are often used 
as tracers of molecular hydrogen. It is valid for molecular clouds in the Galaxy and also 
for external galaxies. The J = 3 level is at 33 K with a critical density of 1.5 x 10 4 cm~ 3 , 
compared to the J = 1 and 2 levels of 12 CO are lying respectively at 5.5 and 17 K above 
the ground level, and with respective critical densities 7.4 x 10 2 and 4.9 x 10 3 cm -3 for 
a kinetic temperature T kin = 30 K in the optically thin limit (Jansen 1995). Therefore, 
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the 12 CO(3-2) transition is a better tracer for the warmer and/or denser molecular gas in 
star forming regions. Furthermore, the ratio of 12 CO(3-2) to (1-0) line emission provides 
a more sensitive measure of the gas temperature and density than the ratio of 12 CO(2-l) 
to (1-0) lines. The first detection of the 12 CO(3-2) line in the nearby starburst galaxy was 
that in IC 342 reported by Ho (1987). Many observations of the 12 CO(1-0) and 12 CO(2-l) 
lines have since been made in local galaxies (e.g. Sanders et al. 1991; Braine & Combes 
1992; Braine et al. 1993; Young et al. 1995; Chini et al. 1996, and many others). Most 
observations of the 12 CO(3-2) line cover the central region of nearby objects (e.g. Devereux 
et al. 1994; Mauersberger et al. 1999; Dumke et al. 2001), where the physical conditions of 
molecular gas may be different from those prevailing in molecular clouds in the disk of the 
Galaxy. Extended maps of 12 CO(3-2) line for three nearby normal galaxies were reported by 
Wielebinski et al. (1999). 

The weaknesses of most of the previous studies on CO molecular gas are the small 
numbers of objects or the lack of uniform selection criteria, whereas their strength is the 
production of maps with high linear resolution. In this paper, we address the aforementioned 
weaknesses by presenting the largest statistical survey of 12 CO (1-0) and 12 CO (3-2) point 
measurements from the central regions in a complete subsample of 60 IR luminous sources 
from SLUGS (Dunne et al. 2000). The angular resolutions of the 12 CO(1-0) observations 
at the Nobeyama Radio Observatory (NRO) and the 12 CO(3-2) observations at the James 
Clerk Maxwell Telescope (JCMT) are nearly identical, so that no assumptions need be made 
about differing beam dilution factors. The beamsize (~ 15") is also nearly identical to that 
of the SLUGS survey by Dunne et al. (2000). Together with the 850 fxm fluxes, plus existing 
data on H I (Thomas et al. 2002), radio continuum (Condon et al. 1990), and far-IR (Dunne 
et al. 2000, NED 1 ), we are able to search for a relationship between the degree of excitation 
of the CO in this sample and the star forming properties. 

It is now suspected that the 12 CO(3-2) line emission makes up a substantial fraction 
of the SCUBA flux at 850 /zm for starburst galaxies (e.g. Papadopoulos & Allen 2000; Zhu 
et al. 2003). Our 12 CO(3-2) observations provide an important database for correcting the 
SCUBA sub-mm data in the SLUGS sample, thus permitting a revised characterization of 
the dust, which will be presented in a separate paper. 

Finally, our investigations on the gas and dust properties in starburst regions provide an 
opportunity to investigate the controversial parameter X, the ratio of molecular hydrogen 
mass M(H 2 ) and the 12 CO(1-0) luminosity L C o, based on the LVG method (Goldreich & 



1 The NASA/IPAC Extragalactic Database (NED) is operated by the Jet propulsion Laboratory, Califor- 
nia Institute of Technology under contract with the National Aeronautics and Space Administration. 
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Kwan 1974). We obtain an estimate for this parameter which depends on Zqo-, the abundance 
of CO relative to H 2 . We also use our estimate of this parameter to investigate the molecular 
gas content of our SLUGS subsample. 

We assume go = 0.5 and H = 75 km s~ x Mpc -1 throughout for consistency with previous 
work on these galaxies (e.g. Dunne et al. 2000). 

2. Observations 

2.1. 12 CO(1-0) Observations at NRO 

Two observing runs measuring the 12 CO(1-0) line at 115.27 GHz for 46 SLUGS galax- 
ies were made with NRO 45m telescope between 2001 April 12 - 19 and 2001 November 
28 - December 1, using the SIS receiver S100 (single side band or SSB). The wide band 
acousto-optical spectrometer (AOS) was used as the back end which provided a total usable 
bandwidth of 250 MHz and a frequency resolution of 250 KHz. At 115.27 GHz, these cor- 
respond to a total velocity coverage of 650 km s" 1 and velocity resolution of 0.65 km s _1 . 
The NRO 45m telescope beamsize (FWHM) at this frequency is 14.6". The standard source 
Mars was used for the determination of the main beam efficiency The measurement yields 
Vmb = 0.34 ± 0.04 after averaging the results from the first and second runs. The average 
uncertainty in the calibration of spectral line is ~ 10%. 

Position switching in azimuth was employed with about 20 seconds each for both ON 
and OFF positions, 10 sec to move from ON to OFF or from OFF back to ON, so that 
about one minute was required for each ON-OFF sequence or scan. Since about 30 scans 
were made for each source, the total integration time per pointing is about 30 minutes. 
Telescope pointing was checked every half hour by observing strong SiO maser lines from 
nearby evolved stars with receiver S40 at 43 GHz. The average pointing error was ~ 3.0" 
(r.m.s.) for the first observing run and ~ 2.6" (r.m.s.) for the second run at low wind 
speed. The weather conditions were average during the NRO observing runs and the system 
temperature for good weather conditions was T sys ~ 535 K. However, half of the observing 
time was lost for each run due to poor weather conditions. The total number of sources 
observed at this transition is 51, out of a total sample of 60. 
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2.2. 12 CO(3-2) Observations at JCMT 

Observations of the 12 CO(3-2) line at 345.80 GHz for the selected 60 SLUGS galaxies 
were obtained with the 15m JCMT on 2001 April 23 - 28 and 2002 June 25. We used the dual- 
channel receiver B3 in single-sideband mode, and the Digital Autocorrelation Spectrometer 
(DAS) spectrometer with a bandwidth of 920 MHz (~ 800 km s _1 ) and a frequency resolution 
of 750 KHz (~ 0.65 km s" 1 ). Focus and pointing were monitored frequently by observing 
bright continuum sources and planets. The average pointing error (r.m.s.) was found to be 
~ 3.2". The observations were performed by using beam-switching at a frequency of 1 Hz 
with a beam throw of 180" in azimuth. The JCMT beamsize (FWHM) at this frequency is 
14.4". The main-beam efficiency was determined by averaging four measurements based on 
observations of Mars, yielding a value of f] m b = 0.50 ± 0.05. The spectral line calibration 
was checked each night by measuring 12 CO(3-2) line emission from IRC+10216. The mean 
calibration error is ~ 10%. The weather conditions for these observations were good, and 
the average system temperature was T sys ~ 600 K. 



3. Results 
3.1. The Spectra 

The reduction of the 12 CO(1-0) data was performed using the NRO NEWSTAR package, 
and the 12 CO J=3-2 data were reduced using the JCMT SPECX package. The 12 CO(1-0) line 
spectra for 51 galaxies together with the spectra of 12 CO(3-2) line emission for 59 galaxies 
are presented in Fig. 1. A linear baseline was subtracted from each spectrum for both data 
sets, and the spectra were smoothed to 10.2 km s _1 and 11.2 km s _1 for the NRO and the 
JCMT data respectively. We converted the measured antenna temperature scale T* of both 
spectra to the main-beam temperature T m b scale by using the relation T m b = T*/r] m b. The 
main-beam efficiencies of NRO and JCMT telescopes are given in § 2. The JCMT and NRO 
line profile shapes mostly agree with each other. The ones that show discrepancies may in 
some instances be due to pointing error. 



3.2. CO and FIR Luminosities 

The velocity-integrated intensities of 12 CO(1-0) line emission I w and 12 CO(3-2) line 
emission 7 32 corresponding to the common telescope beamsize of ~ 15" are shown in Table 1. 
The uncertainties for the integrated line intensities include components due to noise, the un- 
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certainty in the beam efficiencies and the calibration of spectral line, combined quadratically. 
The luminosities associated with the beamsize of ~ 15", in units K km s _1 Mpc 2 Q b , were 
calculated using equations 1(a) and 1(b) below and are presented in Table 2. Both equations 
1(a) and 1(b) were derived by transforming the relevant parameters from the rest frame to 
the observer frame. 

I w Q b , (la) 

(lb) 

where D L is the luminosity distance of a galaxy in Mpc, D L = D (1+z), and D is the 
proper distance of the galaxy (see Table 1), the integrated line intensities Iiq and I32 have 
units of K km s" 1 , and is the beam solid angle. The numerical value of Q b is constant 
for all observations and has been set to unity for simplicity. For conversion to absolute 
luminosities, the value of Q b = 1.1 x 10~ 3 sterradians. 

Fig. 2 shows a comparison between the CO luminosities Leo and FIR luminosities 
Lpjfi for two different FIR color regimes. The quantity Lpj R here corresponds to the volume 
marked by a 15" beam obtained by applying a scale factor to the total FIR luminosity (from 
Dunne et al. 2000). This factor is derived from the original SCUBA images which have a 
resolution of 15" and from radio continuum images (Condon et al. 1990) convolved to 15". 
It measures the ratio of flux within the beam to the total flux. It is less than unity if the 
source is resolved by the telescope beam and equal to unity for an unresolved source. It is 
assumed here that the FIR brightness distribution is similar to that for the sub- mm/cm and 
radio continuum. A linear fit to each transition in the log-log plane yields 

log w Li2 CO{1 _ 0) = (-0.9 ±0.7) + (0.58 ± 0.07)log 10 L FIR , (2a) 
log l0 Li2 CO{3 _ 2) = (-2.4 ±0.5) + (0.70 ±0M)log 10 L FIR (2b) 

The r.m.s. dispersion for the fits is about a factor of 3 and 2 for 12 CO(1-0) and 12 CO(3- 
2) respectively. The scatter is much larger than the errors measured for the telescopes 
indicated by the error bars in the plots. Variations in the CO-to-H 2 conversion factor X 
may be responsible for part of the scatter. The larger scatter seen in the 12 CO(1-0) plot may 
be due to the larger baseline ripple measured at this transition. 
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It is evident that the CO luminosities at both transitions increase non-linearly with 
increasing L FIR , such that the ratio L co j 'L FIR decreases with increasing L FIR . Fig. 3 shows 
a plot of the ratio of Lco/Lfir versus the projected beamsize on the galaxies, to test whether 
this effect is introduced by partially resolving the nearby galaxies. No relationship is evident, 
indicating that the total luminosity of the galaxy is the primary factor. The implications of 
this result will be discussed in § 4.1. 



3.3. The 12 CO(3-2)/(l-0) Line Intensity Ratio 

For the 48 galaxies with data at both transitions, we calculate the intensity ratios of 
the two transitions using the equation, 

?"31 = -W-flO ( 3 ) 

The results along with the asymmetric uncertainties are given in Table 1. The asymmet- 
ric errors of r 3 i are derived analytically, assuming the error distributions of iio and I32 are 
Gaussian functions. The asymmetric confidence intervals in Table 1 represent an enclosed 
probability of 68%. 

The values of r 31 range from 0.22 (for NGC 5665) to 1.72 (for NGC 5937). There are 
seven sources with r 31 > 1, namely UGC 5376, IR 1222-06, UGC 8739, NGC 5433, NGC 
5713, NGC 5937, and MCG +01-42, signifying very high excitation coupled with low optical 
depth. However, UGC 5376, IR 1222-06, and MCG +01-42 have high uncertainties in r 31 , 
leaving only four whose r 31 values are significantly greater than 1.0. In our sample, thirteen 
galaxies have r 3 i ratios less than 0.4, the average value of r 3 i derived from the giant molecular 
clouds in the Galaxy (Sanders et al. 1993). 

The mean value of r 31 for SLUGS sample is 0.66, close to that of 0.63 measured for 28 
galaxies with strong CO emission published by Mauersberger et al. (1999). The quantity 
Lfir ranges from 10 9 to 10 12 L Q in the Mauersberger sample, compared to the range 10 10 - 
10 12 L Q in SLUGS sample. Fig. 4 shows a comparison between the distributions of r 3 i for 
the SLUGS sample and those measured by Mauersberger et al. (1999). The probability of a 
significant difference between these two distributions is ~ 68% using a Kolmogorov-Smirnov 
test. This implies that these two distributions are not significantly different. The galaxies 
observed by Mauersberger et al. (1999) are located between 0.8 and 73 Mpc, much closer, 
and on average less luminous than our objects in the range 25 - 277 Mpc. Thus there is 
no evidence that r 31 depends on the distance of the galaxies over the range 0.8 - 277 Mpc, 
because of the different luminosity ranges represented by these two samples. 
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A search was conducted for relationships between r 3i and a number of parameters 
relating primarily to star formation. The existence of such correlations would be expected 
to shed light on the factors affecting star formation, or the effects of star formation on gas 
properties. These are considered below: 

3.3.1. Projected Beamsize 

Radio continuum maps (Condon et al. 1990), the H I (Thomas et al. 2002), and SCUBA 
850yum maps (Thomas et al. 2002) show that much of the emission in the SLUGS objects 
is extended with respect to the beamsize of ~ 15", especially for the nearer objects. On 
average we detect about 45% of all the emission from these galaxies assuming the CO is 
distributed in a similar way to the 850 /xm and radio continuum emission. Therefore, the 
CO measurements may reflect the physical conditions in the center regions for the galaxies 
at lower cz. Maps of 12 CO(3-2) line emission in 12 nearby galaxies by Dumke et al. (2001) 
have shown that the 12 CO(3-2)/(l-0) line intensity ratio decreases outward from the centers 
of many galaxies. In principle, for a given beamsize, the average line ratios obtained from 
the single pixel measurement should therefore decrease with increasing distance, since more 
emission from lesser excited CO lines arising from the disk regions is detected by the telescope 
beam. Fig. 5 shows a plot of r 31 versus projected beam size. Similar to the relation between 
the ratio of Lco/Lfir and the linear size sampled (see Fig. 3), no correlation between these 
quantities is evident. The lack of such a correlation may be due to sample selection effect, 
however, since the more distant galaxies may have a larger fraction of highly excited CO 
molecular gas. 

3.3.2. Dust and Gas 

Large amounts of FIR emission are generated by the dust heated by young massive 
stars in starburst systems, in which case the FIR luminosity is a good measurement of the 
star formation rate (SFR) (Kennicutt 1998a). Fig. 6 shows the 12 CO line ratio versus dust 
temperature T dust and FIR luminosity L FIR within a 15" beam. The dust temperatures 
are from Dunne et al. (2000) and the FIR luminosities are scaled to a 15" beam using 
sub-mm and radio maps as described in § 3.2. No significant correlation is evident in these 
plots. Other SFR indicators are radio continuum luminosity L ra( ao and the optical color 
index (B — V). Plots for these are not shown, but no obvious correlation was found with 
the parameter r 31 . In Fig. 7, we examine the dependence of r 31 on dust mass M^ust and 
molecular hydrogen mass M(H 2 ). The dust masses were obtained from Dunne et al. (2000) 
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and scaled to a 15" beam as described earlier. The molecular gas masses for the SLUGS 
sample are derived from our CO luminosities by applying the conversion factor X = 2.7 x 
10 19 cm" 2 [K km s" 1 ]" 1 derived in § 4.3.1. The result is M(H 2 ) = 1.1 x 10 3 D 2 L /(l+z) S co 
M Q by scaling the result of Kenney & Young (1989), where Dl is the luminosity distance 
of a galaxy in Mpc, and Sco is the 12 CO(1-0) flux measured within a 15" beam, Sco = 2.1 
iio Jy km s _1 . The values of M(H 2 ) are also presented in Table 2. Again, no relationship is 
evident. 



3.3.3. Surface Brightnesses and Color Indices 

The average surface brightness for all wavebands is defined here as the ratio of the 
total flux to the square of the optical diameter. Thus it measures the surface brightness 
averaged over the optical diameter, and provides a measure of the concentration of star 
forming activity. Fig. 8 shows a plot of r 31 versus the surface brightness for mid-IR, far- 
IR, sub-mm, and H I emission. We find two possible correlations, one is between r 31 and 
FIR surface brightness, and another one is between the r 3 i and SCUBA 850 /im surface 



brightness. Fig. 9 shows the line ratio versus different color indices or ratios: 
No significant correlation is found between these ratios. 
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3.3.4- Star Formation Efficiency 

The ratio Lfir/M(H.2) is often used as an indicator of the star formation efficiency 
(SFE) which gives the star formation rate per unit mass of molecular gas. In Table 2, we 
present the results of SFE calculations for our SLUGS sample. Here the mass of molecular 
gas has been derived in the same manner as that in § 3.3.2. These values are an order 
of magnitude higher than what would be obtained if the conventional value of X is used. 
Fig. 10 shows a significant correlation between r 3i and the SFE within the telescope beam. 
The sample shown in Fig. 10 is divided into two ranges by gas mass centered at M(H 2 ) = 
10 8 M , and also by dust FIR luminosity centered at Lfir = 10 10 L Q which are indicated 
by four different symbols in the plot. The linear correlation coefficient is obtained from the 
data with SFE < 200 L Q /M . The coefficient is 0.74 at a significance level of 1.4 x 10~ 7 
(i.e. probability that r 31 and SFE are unrelated). However, the correlation is diminished at 
SFE > 200 L Q /M , where the ratios spread between 0.5 and 1.72. The significance of the 
dependence of r 31 on SFE will be discussed in § 4.2. 
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4. 



Discussion 



4.1. 



The Relation Between CO and FIR Luminosities 



As noted in § 3.2, the non-linear variation of CO with FIR luminosity implies a decrease 
in the ratio L C o/ L FIR with increasing L FIR (see Fig. 2). For 12 CO(1-0) this ratio decreases 
by more than an order of magnitude over the range L FIR = 5 x 10 9 - 10 12 L . Part of this 
variation might be attributable to the effect of beamsize if the ratio is dependent on the 
fraction of the galaxy resolved by the beam. However, as noted in § 3.2, Fig. 3 shows that 
the beam effects do not seem to influence the Lco/Lfir ratio, and that for any plausible 
gas excitation gradient (i.e. hotter gas confined in the center) such effects would show in 
the r 31 ratio as well, but Fig. 5 shows no evidence that more distant objects have smaller 
values of r 31 (see § 3.3.1). In addition, this effect would be expected to produce an increase 
rather than a decrease in Lco/Lfir with increasing Lfir because of the strong sensitivity 
of Lfir to dust temperature. Therefore this ratio depends intrinsically on the total FIR 
luminosity. It should be noted that our results imply that for high z galaxies where L FIR is 
exceedingly high, the corresponding values of L C o will be seriously overestimated if a linear 
relation between Leo and Lfir were assumed. 

A possible cause of the variation in the Lco/Lfir ratio is an increase in Lfir per unit 
mass of dust with increasing Lfir- To investigate this point further, we first plot in Fig. 11 
L C o versus M dust to determine this relationship is linear. Linear fits to these relations in the 
log-log plane yield 



We see indeed that the slopes are closer to unity than with logio-^co versus Xogi^LpiR 
(Equation (2a) and (2b)), though there is still a significant departures from a linear relation 
for I/i2(7O(i-0)- Secondly, we note that Fig. 2 shows a trend that the Lco/Lfir ratio decreases 
with increasing dust temperature which is indicated by the two different FIR color ranges. 
This implies that the galaxies of the SLUGS sample which have higher luminosity is because 
they have a higher temperature in addition to simply more dust and gas, which would alone 
yield a linear relation between L C o an d L F ir- The correlation found between L 60 ^ m and dust 
temperature by Dunne et al. (2000) also support this argument. This effect is consistent with 
the systematically higher SFE's deduced for such objects, with the warmer IRAS colors, 
and with the saturation of the r 31 ratio at high SFE's (see Fig. 10). 



log w Li2 C o(i-o) 
log w Li2 C o(z-2) 



(-0.4 ± 0.6) + (0.77 ± 0.09)log w M dust 
(-1.8 ± 0.5) + (0.95 ± 0.07)log w M dust 



(4a) 
(4b) 
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As noted in § 3.3.2, L FIR is widely used as an indicator of the SFR, and M(H 2 ) is 
often computed from the CO flux. Thus, the non-linear relation between L C o and L FIR 
may reflect the behavior of star formation law for the nearby IR luminous galaxies. The 
star formation law can be represented as a power law that can be expressed in terms of 
the star formation rate per unit area Y^sfr and the gas density S 9as , i.e. Hsfr = 
first introduced by Schmidt (1959). Previous studies have shown that the activity of disk- 
averaged star formation is correlated with the mean total molecular gas content (Kennicutt 
1989; Boselli 1994; Boselli et al. 1995), and this correlation is found to be stronger with H 2 
gas in starburst regions (Kennicutt 1998b). Fig. 12 shows a possible correlation between the 
SFR per unit area and the H 2 gas density £# 2 measured within a 15" beam for our SLUGS 
sample. Both T>sfr and E# 2 are corrected for the inclination of galaxies in the SLUGS 
sample. The correction for the galaxy inclination is necessary because the Schmidt Law 
involves T,$fr and E# 2 when viewed normal to the disk. The disperson in the fit is about a 
factor of 3. The slope based on a linear fit in the log-log plane yields 1.4 ± 0.3. This value is 
consistent with the Schmidt law applied to the circumnuclear starbursts (Kennicutt 1998b), 
though the uncertainty in the slope is too high to obtain a definitive value for the index. 
A possible reason for this is that our result refers to different regions in different galaxies, 
depending on the projected beamsize, whereas the work by Kennicutt and others consistently 
refer to disk- averaged parameters. The Schmidt Law involves H I+H 2 , but ignoring H I gas, 
which averages about 30% of the total gas mass measured within a 15" beam (see § 4.3.2), 
is most likely not the cause of the scatter seen in Fig. 12, since a study by Wong & Blitz 
(2002) with a high angular resolution shows no correlation between H I and star formation 
rate. 



4.2. Correlation between r 31 and Star Formation Parameters 

We examined in § 3.3 the excitation of CO and its relation with the properties of gas/dust 
and star formation in the central starburst regions in SLUGS sample. There are no significant 
correlations between r 31 and the distance of galaxies, star formation rate, dust temperature 
and mass, H 2 gas mass, the color indices, and the luminosities of IR and radio continuum. 
The lack of correlation of r 31 with properties relating to total star formation reflects a range 
of localized conditions in the molecular clouds. The possible dependence of r 3i on FIR/sub- 
mm surface brightness (see Fig. 8) and SFE (see Fig. 10) suggests that the higher degree 
of CO excitation is related to a higher concentration and efficiency of star forming activity. 
Such conditions would arise in an intense starburst where the surface density of such activity 
is high, consuming most of the gas in the region. The saturation effect of r 3 i seen at higher 
SFE which is probably due to thermalization of the CO levels, reflecting denser and warmer 
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gas in regions of high star formation efficiency. This interpretation is consistent also with the 
non-linear relation between L C o an d L FIR (see Fig. 2), and with higher color temperatures 
associated with the most IR luminous and most efficient star forming galaxies. It is not clear 
whether these correlations reflect a high excitation ratio as the underlying cause or the effect 
of star forming activity. 

The regime where r 31 > 1.0 requires optically thin CO lines from extremely warm and/or 
dense gas. The relatively small number of galaxies with high values of r 31 might be associated 
with the highly excited gas from a short-lived gas phase in the starburst systems. The lower 
line ratios may be associated with gas that is cooler and less denser. 

4.3. The Molecular Gas Content 

The investigation of the gas and dust content depends critically on a reliable estimate of 
the controversial parameter X. Studies have shown that this parameter varies from galaxy 
to galaxy (Booth & Aalto 1998; Boselli et al. 2002), and it is thought to be higher in metal- 
poor galaxies and lower in starburst galaxies than in Galactic molecular clouds, where X 
is about 2.8 x 10 20 cm" 2 [K km s" 1 ]" 1 (Bloemen et al. 1986; Strong et al. 1988). Thus in 
starburst galaxies, application of the standard factor can produce a significant overestimate 
of molecular hydrogen mass (Solomon et al. 1997; Downes & Solomon 1998; Zhu et al. 2003). 

4.3.1. The CO-to-H 2 Conversion Factor X from An LVG Method 

In this section, we examine the molecular gas content and properties using large velocity 
gradient (LVG) models. The LVG method is capable in principle of yielding the column 
density of CO, and hence of H 2 , if an abundance ratio Z C o — [CO/H 2 ] is assumed. Hence, 
this method should yield an estimate of the X factor independent of other considerations 
such as gas-to-dust ratio, which is discussed in the next section. We assume that the ISM 
comprises molecular clouds characterized by a single set of average physical properties. It 
is well understood that models incorporating more data for starburst galaxies require at 
least two components - a prevailing optically thin warm diffuse component, and an optically 
thick cool dense component with a smaller filling factor (e.g. Aalto et al. 1995). However, 
comparisons between single and double component models yield approximately the same 
column density for H 2 , which is weighted strongly toward the prevailing diffuse component 
(e.g. Zhu et al. 2003). 

The LVG approach can be formulated into an equation yielding the X factor from the 
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LVG parameters given by 

X = f^- (5, 

^CO^rad 

where n(H 2 ) is the H 2 gas density, T rad is the radiation temperature for the 12 CO(1-0) 
line transition, A = ^ v C /dr) ' an< ^ (dv/dr) is the velocity gradient associated with the individual 
clouds (see Appendix A). There are not enough data to uniquely determine the parameters 
in this equation, which usually require multiple line excitation ratios, measured brightness 
temperatures, and isotope intensity ratios to fully constrain the conditions in the molecular 
gas. In our case, it is therefore necessary to make some reasonable assumptions to constrain 
the LVG solutions. 

The highest degree of uncertainty concerns the parameter Zco which is unknown and 
needs to be addressed. Values widely quoted for starburst galaxies are in the range 10~ 5 - 
10~ 4 (e.g. Booth & Aalto 1998; Mao et al. 2000). Unfortunately, there are few if any reliable 
determinations for starburst galaxies. Within the Milky Way, measurements of Zco an d 
chemical models for dark and star forming clouds yield values within the range 5 x 10~ 5 
- 2.7 x 10" 4 (Blake et al. 1987; Farquhar et al. 1997; Hartquist et al. 1998; van Dishoeck 
1998; van Dishoeck & Blake 1998). It seems unlikely therefore that Z C o < 10~ 5 for starburst 
galaxies unless their metallicity is unusually low, which would occur only in low luminosity 
systems. For the SLUGS sample, therefore, we adopt Zco = 10~ 4 , bearing in mind the 
dependence of our computed values for X on this parameter. 

We adopt dust temperatures for each galaxy from Dunne et al. (2000), assuming T dust 
to be representative of gas kinetic temperature T kin . Fig. 13 shows a series of plots of r 3i 
versus T dust for our subsample combined with theoretical curves based on LVG models. The 
theoretical curves represent r 31 versus T kin for a variety of assumed densities n(H 2 ) based 
on a series of LVG models using Zco — 10~ 4 , and values of A = 10~ 7 , 10" 6 , 10~ 5 , and 10~ 4 
(km s" 1 pc" 1 )" 1 . In each case it is assumed that T^n = T dust . Note the convergence of the 
curves at higher densities reflecting the saturation of the line ratios due to thermalization of 
the optically thick CO transitions. Note also that the ratios r 31 > 1.0 can not be readily fit 
by the models with A > 10~ 6 (km s" 1 pc" 1 )" 1 , a point to which we return later. 

Thus for each galaxy in one of the data plots, it is possible to derive a molecular gas 
density and radiation temperature T rad which is dependent upon A. Using Equation (5), 
we then compute a range of X parameter for the IR luminous galaxies. For each A, only 
points that fall below the upper limits of r 31 due to saturation are used for computing the 
distribution of X. For example, for A = 10~ 5 (km s" 1 pc" 1 )" 1 , thermalization occurs at r 3i = 
0.86, corresponding to T dust = 55 K, n(H 2 ) = 2.4 x 10 3 cm" 3 , and T rad = 27 K. The resulting 
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range of X is 1 - 3 x 10 19 cm -2 [K km s" 1 ]" 1 . For points falling above the saturation curves 
(r 31 > 0.86) for A = 10~ 5 (km s" 1 pc" 1 )" 1 (see Fig. 13), we obtain X values in a range 2.4 x 
10 19 - 6.7 x 10 19 cm" 2 [K km s" 1 ]" 1 from the plots corresponding to A = 10" 6 and 10~ 7 (km 
s" 1 pc" 1 )" 1 . However, the latter value of A is virtually unacceptable, since it corresponds 
to an impossibly high velocity gradient or impossibly low CO abundance. Therefore, Tfc; n 
would have to be higher than the measured Td us t for those points to lie below the saturation 
curves for A < 1CT 6 (km s" 1 pc" 1 )" 1 . Fig. 14 shows the resulting distribution of the derived 
X and its average values for different A. It is clear that regardless of A, the values of X 
derived from starburst galaxies are systematically lower than the value X = 2.8 x 10 20 cm" 2 
[K km s" 1 ]" 1 derived from the disk of the Galaxy. 

In order to further constrain the range of the X factor and A shown in Fig. 14, we 
compare the CO isotope intensity ratios R 1Q = J[ 12 CO(1-0)]//[ 13 CO(1-0)] predicted by the 



LVG model with observed values. We adopted the isotope abundance ratio as 



2 CO 



40 



13 CO 

75 (Wilson & Rood 1994). The CO isotope intensity ratios data are taken from Aalto et al. 
(1995) and Taniguchi et al. (1999) for 34 IR luminous galaxies which have a range of FIR 
luminosities similar to the SLUGS sample, i.e. L FIR > 10 10 L©. Fig. 15 shows a comparison 
between the distributions of computed and observed isotope ratios for the sample as a whole 
for each value of A, and corresponding value for the mean derived X factor. The plots show 
that the best agreement with the observed distribution of Riq is obtained for A = 10~ 5 (km 
s" 1 pc" 1 )" 1 , corresponding to an average value for X = 2.7 x 10 19 cm" 2 [K km s" 1 ]" 1 . The 
other plots assuming other values of A show marked disagreement with the observed isotope 
line ratios. 

The resulting value for X (2.7 x 10 19 cm" 2 [K km s" 1 ]" 1 ) in starburst regions is thus 
ten times smaller than the conventional value of X for GMC's in the Milky Way. But the 
result is comparable with X = 0.5 x 10 20 cm" 2 [K km s" 1 ]" 1 estimated from diffuse clouds 
in the Galaxy by Polk et al. (1988), and with that found for extreme starbursts in nearby 
galaxies by Downes & Solomon (1998) and for starburst galaxy M 82 by Mao et al. (2000). 

We repeated the LVG model calculations for a higher CO isotope abundance ratio of 75 
(Wilson & Rood 1994), and found that the average value of X for Z C o — 10~ 4 becomes 8.2 
x 10 19 cm -2 [K km s" 1 ]" 1 . We also repeated the entire analysis by increasing the adopted 
kinetic temperatures by 50%, and found that the resulting value of X remains essentially the 
same, since for a given A, the lower densities required are largely offset by lower values of the 
radiation temperature (see Equation (5)). Our LVG analysis also shows that the predicted 
isotope ratios 12 CO(2-l)/ 13 CO(2-l) are in satisfactory agreement with the observed ratios 
reported by Taniguchi et al. (1999). The results do not seem to support the conclusion by 
Taniguchi et al. (1999) that the comparison between CO(2-l) and CO(1-0) isotope intensity 
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ratios may imply that the unusually large values in starburst systems indicate an isotope 
abundance anomaly in starburst galaxies. Our results are consistent with effects of excitation 
and optical depth. 

4.3.2. TheM (H 2 )/M (H I) Ratio 

In this section, we briefly consider a small subset of our sample with H I maps available 
to assess the combined molecular and atomic gas content within the telescope beam. The 
subset comprises twelve galaxies with an average L FIR = 4.1 x 10 10 L Q (Thomas et al. 2002). 
The H I masses for our 15" beam were estimated from the peak H I column densities maps 
at 25" resolution by assuming a uniform H I brightness distribution within the 25" beam. 
The mean (± r.m.s.) projected beamsize for the subset is 3.6 (± 1.1) kpc. 

The overall mean (± r.m.s.) values for M(H 2 )/M(H I) ratio and combined H 2 +H I 
face-on surface density are respectively 2.0 ±1.5 and 42 ± 18 M pc~ 2 . These quantities 
must be considered respectively as upper and lower limits however, since the H I masses may 
be underestimated due to the effects of absorption by H I against the radio continuum in 
the nuclear regions (Thomas et al. 2002). This is particularly the case for one galaxy, NGC 
5900, for which M(H 2 )/M(H I) = 21, and this object was omitted from the mean value 
quoted above. The mean ratio M(H 2 )/M(H I) is substantially larger than that found for 
optically selected late type galaxies (0.15) by Boselli et al (2002). This is partly attributable 
to the fact that their values are global measures, whereas our values refer to the nuclear 
regions of galaxies selected according to high star formation rates. The H 2 +H I surface 
density is intermediate between the central surface densities for galaxies with lower star 
forming rates like the Milky Way, which is about 10 M pc -2 (Elmegreen 1989) and those 
for ultra-luminous IR galaxies, which are up to 10 5 M Q pc~ 2 (Aalto et al. 1995). 

It should be emphasized that these estimates refer to comparatively warm gas heated 
by star formation in these systems. The outer disks of galaxies may be hiding substan- 
tial amounts of cool H 2 , where the radiative efficiency of CO is low (i.e. X is high). The 
M(H.2)/Md U st ratio will be discussed in a separate paper (in preparation), where the dust 
mass of the SLUGS sample (Dunne et al. 2000) will be corrected for the effects of contami- 
nation of the SUCBA 850 /iin filter. 



-16- 



4-3.3. Virial Stability of The Molecular Clouds 

Using the virial theorem, the velocity gradient under the virial condition for a spherical 
cloud with a mean gas density < n > is given by, 



(dv\ 
\dr) 



VIR 



0.65a 2 f— 5 — ^—\ 2 kms 1 pc 1 (6) 



where a is in the range 0.5 - 2.5 (Papadopoulos 1998; Bryant & Scoville 1996) depending 
primarily on the assumed density profile. For a mean gas density of 1.4 x 10 3 cm -3 for r 31 
< 0.8 derived from the LVG model for our IR luminous galaxies (see Fig. 13) with A = 
10 -5 (km s" 1 pc -1 ) -1 and Z C o — 10~ 4 , we estimate a virial velocity gradient of 1.1 km s~ x 
pc~ 4 for a = 2.5. This virial value is about 9-90 times smaller than the velocity gradient 
of 10 - 100 km s~ 4 pc _1 also derived from LVG modeling using A = 10~ 5 - 10~ 6 (km s _1 
pc -1 ) -1 and Zco = 10~ 4 . The same calculation was made for a lower r 3 i using the source 
UGC 6436 with r 31 = 0.35, and a corresponding gas density of 7.5 x 10 2 cm~ 3 . It was 

found that is still about 11 times smaller than the LVG velocity gradient. Since 

\ dr J VIR 

such low values of r 3 i are more sensitive to density than more saturated values near unity, 
this result provides a strong confirmation of the non-virial expansion velocities. As noted 
earlier, there are a number of galaxies with high 12 CO line ratios (r 31 ~ 0.8 - 1.72) seen 
in Fig. 13, where A may be lower than 10 -5 (km s~ 4 pc -1 ) -1 . The above estimates imply 
that the molecular clouds in the starburst regions are not gravitationally bound, confirming 
the suggestion recently made by Zhu et al. (2003, ApJ in press) for the Antennae galaxies, 
unless one is willing to accept a significantly lower value for Z C o- The cause could be the 
destruction by the stellar winds and expanding shells of new-born massive stars. We note 
that the non-virialized component emanates from the gas phase that dominates the 12 CO 
emission, and not necessarily from the entire cloud mass, in accordance with two-phase gas 
models (Aalto et al. 1995; Papadopoulos & Seaquist 1999). 



5. Conclusions 



The main conclusions of this work are: 

1. The Lco/Lfir ratios measured from the SLUGS sample decrease with increasing 
FIR luminosity. The non-linearity between L C o/Lfir and L FIR implies that the more IR 
luminous galaxies have higher dust temperatures than those at the low luminosity end. This 
non-linear relation was found to be consistent with that expect for the Schmidt Law for star 
formation with an index N = 1.4 ± 0.3, though the index is too uncertain to provide a new 
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and definitive Schmidt relation. 

2. The degree of CO excitation measured by the 12 CO(3-2)/(l-0) line intensity ratios 
vary greatly from galaxy to galaxy in the SLUGS sample. There is a trend for the r 31 
ratios to increase with increasing concentration and efficiency of star forming activity. The 
saturation of r 3i seen at higher SFE implies that the gas is denser and warmer in regions 
of high SFE which is consistent with the non-linearity between L C o/Lfir and L FIR . 

3. Our 12 CO line measurements together with dust and CO isotope data taken from 
the literature are modeled using the LVG method to estimate that the CO-to-H 2 conversion 
factor lies between 2.7 x 10 19 and 6.6 x 10 19 cm~ 2 [K km s -1 ]" 1 for SLUGS galaxies, which 
is about 4-10 times lower than the conventional X derived from the Galaxy. 

4. For a subset of 12 galaxies with available H I maps, the total gas face-on surface 
density within about 2 kpc of the center averages about 42 M pc~ 2 , which is intermediate 
between that for disk galaxies with low and high star formation rates. 

5. Our Virial analysis shows that the molecular clouds in starburst regions are not 
gravitationally bound unless one is willing to accept a 9 - 90 times lower [CO/H 2 ] abundance 
ratio. Most of the 12 CO line emission originates from the non-virialized warm and diffuse 
gas clouds. 
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A. Derivation of The CO-to-H 2 Conversion Factor X From LVG Model 

From the LVG model for a spherical molecular cloud, the column density per unit 
velocity range in the average cloud is given by 



(Al) 
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where n(H 2 ) is the density of H 2 , A = J^;, Z C o is the CO abundance ratio [CO/H 2 ], 
and dv/dr is the velocity gradient within the cloud. The total CO column density per cloud 
Nco, integrated over velocity is then 

N co = n{H 2 )AAv (A2) 
where Av is the velocity width of an individual cloud. 



The column density N C o averaged over the beam area can be written as 



or, 



N co = n(H 2 ) AAv—^—— (A3) 



N doud a Av 



N co = n(H 2 )A -P^tttt AV (A4) 

L ^beam u ' J 

where N doud is the number of clouds in the beam, a is the cloud cross section area, 
^■beam is the beam area, and Ay is the line width of all emission detected in the beam. 

The quantity inside inside the brackets of Equation (A4) is the brightness dilution factor 
representing the product of the geometric and velocity filling factors within the beam. 

Assuming for simplicity, a rectangular line profile, then 



Nco = n(H 2 )A AV (A5) 

- rad J 

where the radiation temperature T rad is calculated from LVG modeling for 12 CO (1-0) 
or (3-2) line emission, and T m b is the corresponding main-beam temperature. 

The quantity T mb AV can be replaced by J T mb dV, the line flux in units of brightness. 
Thus 



% = ^/v (A6) 

J- rad J 

Setting N C o — Z C o N(H 2 ), where N(H 2 ) is the corresponding column density of H 2 , 
and X = N( r H ^ ; where 
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ho = J T mb dV (A7) 
for 12 CO (1-0) line emission. We obtain 

X - (AS) 
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Fig. 1. — Spectra of 12 CO (1-0) and 12 CO (3-2) for individual sources. The velocity scale is 
Heliocentric. The temperature scale is in units of main beam temperature. The letter p 
denotes those galaxies that belong to a pair. 
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Fig. 1. — (continued) 
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Fig. 1. — (continued) 
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Fig. 1. — (continued) 
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Fig. 2. — Plots of CO luminosity versus FIR luminosity, both associated with the 15" beam 
for (a) 12 CO(1-0) and (b) 12 CO(3-2) emission. The lines represent linear regression fits to 
the data. The data points have been segregated according to two FIR color regimes in the 
plot. See text for details. 
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Fig. 3. — Plots of the ratio of CO-to-FIR luminosity versus the projected beamsize for 
12 CO(1-0) and (b) 12 CO(3-2) emission. 
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Fig. 4. — Distribution of the ratio r$\ for IR luminous galaxies, (a) by Mauersberger et al. 
(1999); (b) this paper. 
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Fig. 5. — Plot of r 3 i versus the projected beamsize for the SLUGS subsample. 
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Fig. 6. — Plots of (a) r 31 versus the dust temperature T dust , (b) r 31 versus the FIR luminosity 
Lfir applicable to the 15" beam. 
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Fig. 7. — Plots of (a) r 31 versus dust mass M dust , (b) r 31 versus molecular gas mass M(H 2 ) 
applicable to the 15" beam. 
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Fig. 8. — Plots showing r 31 versus the surface brightness averaged over the optical diameter 
for (a) mid-IR 14.3 fxm, (b) far-IR 60 fjm, (c) SCUBA 850 /im, and (d) H I 21cm line. 
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Fig. 9. — Plots of r 31 versus the color indices or ratios involving (a) far-IR and (b) sub-mm 
to H I 21cm emission. 
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Fig. 10. — Plot of versus the star formation efficiency SFE = Lfir/M(B.2) measured 
within a 15" beam. The line represents a linear regression fit to the data with SFE < 200 
L /M . The data points have been segregated according to gas mass-class and dust IR 
luminosity-class to show the relationship between L FIR and M(H 2 ) and position in the plot. 
See text for details. 
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Fig. 11. — Plots of CO luminosities for the 15" beam versus dust mass from Dunne et 
al. (2000) for (a) 12 CO(1-0) and (b) 12 CO(3-2) emission. The dust masses have not been 
corrected for contamination of the SCUBA 850 fim filter by 12 CO(3-2) line emission. 
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Fig. 13. — Plots of r 31 versus dust temperature T^ ust for different densities of molecular 
hydrogen logio^H^) based on LVG model computations for four A values between 10~ 7 and 
10~ 4 (km s -1 pc -1 ) -1 . Also shown are the observed data, assuming Tkm = T^ust- The 12 CO 
abundance Zco is fixed at 10~ 4 , and the isotope abundance ratio 
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Fig. 14. — Distribution of the X conversion factor derived from LVG computations and its 
mean value X M ean for different values of A = ■ The units of X is cm _2 [K km s -1 ] -1 

and A is in (km s _1 pc -1 ) -1 . 
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Fig. 15. — Histograms of the 12 CO(1-0)/ 13 CO(1-0) intensity ratio R 10 for different values of 
A. Dotted lines: LVG modeling results from this work. Solid lines: data taken from Aalto 
et al. (1995) and Taniguchi et al. (1999). The ^gg] abundance ratio is fixed at 40. The 
values for R 10 observed correspond to a galaxy sample with L FIR > 10 10 L Q , similar to the 
SLUGS subsample. 



Table 1. 12 C0 line parameters of observed 60 IR luminous SLUGS galaxies. 
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a Galaxy name, the letter p denotes those galaxies that belong to a pair. 
b Right ascension ,72000 epoch from NED 
c Declination J2000 epoch from NED. 
d Recession velocity from NED. 

°Galaxy proper distance to the object at the present time. 

f The 12 CO(1-0) line intensity Iiq — J T^^dv in K km s — 1 . Upper limits refer to 3 times the standard error due to the noise. 
^Standard error of Iiq. 

h Thc 12 CO(3-2) line intensity J32 — / T mb dv in K km s — 1 . Upper limits refer to 3 times the standard error due to the noise. 
1 Standard error of 1^2- 
J Thc line intensity ratio r^i- 
k The negative uncertainty in r^i . 
^he positive uncertainty in r^\. 
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Table 2. CO luminosities and masses. 
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8.78 


27 


NGC 5258p 


6.63 




5.47 




5.10 




8.84 


24 


UGC 8739 


4.90 




5.28 




5.32 




8.65 


73 


NGC 5371 


2.48 








4.25 








NGC 5394p 


3.38 




5.32 




4.87 




8.69 


15 


NGC 5395p 


3.39 




4.81 








8.17 


8 


NGC 5433 


4.24 




4.89 




5.12 




8.26 


160 


NGC 5426p 


2.55 




4.34 




3.77 




7.70 


14 


NGC 5427p 


2.60 




4.18 




3.98 








ZW 247.020 


7.48 




5.45 




5.40 




8.82 


169 


NGC 5600 


2.25 




4.01 




3.83 




7.37 


121 


NGC 5653 


3.46 




5.13 




4.92 




8.49 


68 


NGC 5665 


2.16 




4.83 




4.16 




8.19 


11 


NGC 5676 


2.05 








4.06 








NGC 5713 


1.85 




4.02 




4.24 




7.38 


282 


UGC 9618p 


9.68 








4.92 








NGC 5792 


1.87 




4.36 




4.31 




7.72 


127 


ZW 049.057 


3.79 




4.85 




4.77 




8.22 


517 


NGC 5900 


2.44 




4.91 




4.57 




8.28 


53 


1 ZW 107 


11.79 








5.53 








NGC 5929p 


2.42 




3.80 




3.74 








NGC 5930p 


2.60 




4.43 




4.16 




7.79 


99 



-44- 



Table 2 — Continued 



Source a 


P.B. b 


logio 


Ll2 CO ( 1 _ 


C 

0) 


log io 


il2 CO(3-2) d 


logi M(H 2 ) 


SFE { 




(kpc) 


(K km 


s" 1 Mpc 2 


n 6 ) 


(K km 


s- 1 Mpc 2 n b ) 


(M ) 


(Lq/Mq) 


IR 1525+36 


16.31 










5.64 






NGC 5936 


3.90 




5.19 






5.02 


8.56 


82 


NGC 5937 


2.73 




4.41 






4.64 


7.77 


469 


NGC 5953p 


1.91 




4.49 






4.42 


7.85 


58 


NGC 5954p 


1.90 




4.41 






3.99 


7.78 


61 


ARP 220 


5.29 




5.80 






5.57 


9.17 


467 


IR 1533-05 


7.99 










5.00 






NGC 5962 


1.90 




4.76 






4.26 


8.12 


13 


NGC 5990 


3.73 




5.34 






4.85 


8.71 


58 


NGC 6052 


4.59 




5.07 






4.83 


8.44 


102 


MCG +01-42 


7.16 




5.38 






5.40 


8.75 


247 


NGC 6181 


2.31 




4.77 






4.32 


8.13 


26 



a Galaxy name, the letter p denotes those galaxies that belong to a pair. 
b The projected beamsize (P.B.) at the source. 

c The 12 CO(1-0) luminosity measured within a 15" beam, Q b = 1 beam area. 
d The 12 CO(3-2) luminosity measured within a 15" beam, fl b = 1 beam area. 
e Molecular hydrogen gas mass measured within a 15" beam. 
f Star formation efficiency measured within a 15" beam. 



